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ABSTRACT 

We present observations of Ha emission line profiles taken at Magellan Obser- 
vatory for a sample of 39 young low-mass stars and brown dwarfs in the Taurus 
and Chamaeleon I star forming regions. We have identified 11 new substellar ac- 
cretors, more than tripling the number of known brown dwarfs with measurable 
accretion activity. These include the lowest-mass objects yet seen with accretion, 
with masses down to ~ 0.015 Mq. Using models of Ha emission produced in 
magnetospheric accretion flows, the most widely applicable primary calibrator 
now available, we determine the first estimates of mass accretion rates for ob- 
jects at such extremely low masses. For the six objects with masses < 0.03 Mq, 
we find accretion rates of ~ 5 x 10^^^ MqJt~^, among the smallest yet mea- 
sured. These new results continue the trend of decreasing mass accretion rate 
with decreasing (sub) stellar mass that we have noted previously for samples of 
more massive objects; the overall correlation is M oc M^'^, and now extends over 
a mass range of over two orders of magnitude. Finally, the absence of a disconti- 
nuity in the distribution of accretion rates with mass tends to suggest that stars 
and brown dwarfs share similar formation histories. 

Subject headings: accretion disks, brown dwarfs, stars: emission-line, pre-main 
sequence, circumstellar matter. 
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1. Introduction 

The subject of the formation of substellar objects has garnered significant attention in 
the past few years. It remains unclear whether brown dwarfs form from the collapse of very 
low mass cloud cores, such as produced by turbulent fragmentation (Padoan & Nordlund 
2002), or whether they are created via ejection of low- mass embryos from the main reservoirs 
of gas in a molecular cloud (Reipurth & Clarke 2001; Bate et al. 2002). These two possibilities 
can be investigated by comparing the characteristics of newborn substellar objects to those 
of their T Tauri counterparts at stellar masses. 

The paradigm for low-mass star formation begins with the collapse of a molecular cloud 
core onto a protostar and circumstellar accretion disk (Adams et al. 1987). The protostar 
continues to grow via magnetospheric accretion, whereby the inner regions of the circumstel- 
lar disk are disrupted by the stellar magnetosphere, which then channels viscously accreting 
material out of the disk plane and towards the star along magnetic field lines (Konigl 1991; 
Shu et al. 1994). Signposts of this process, which have been directly observed in classical 
T Tauri stars (hereafter, CTTSs), include: infrared (IR) emission from dust at a range of 
temperatures in the disk heated by stellar irradiation and viscous dissipation (e.g., Meyer 
et al. (1997); MuzeroUe et al. (2003a)); blue/UV continuum excess emission from the accre- 
tion shock formed as accreting material falls onto the stellar surface (Hartigan et al. 1995; 
Valenti et al. 1993; GuUbring et al. 1998; Calvet et al. 1998); broad permitted emission fines 
produced in the ballistic magnetospheric gas flows (Muzerolle et al. (2001), hereafter MCH); 
forbidden emission lines produced in accretion-powered winds and jets (Hartigan et al. 1995). 

Recent studies have begun to examine these various signposts near and below the hy- 
drogen burning limit to ascertain whether the same formation mechanism applies at both 
stellar and substellar masses. High-resolution optical spectroscopy of young brown dwarfs has 
indeed provided considerable evidence for magnetically-mediated disk accretion, including 
broad permitted emission line profiles (especially Ha ) and optical continuum veiling (White 
& Basri (2003), hereafter, WB03; Muzerolle et al. (2003b), hereafter, M03; Jayawardhana 
et al. (2003b)). Furthermore, detection of IR excess emission in many objects confirms the 
presence of irradiated circumstellar disks (Comeron et al. 1998, 2000; Luhman 1999; Muench 
et al. 2001; Natta et al. 2002; Jayawardhana et al. 2003a). 

In this paper, we present high-resolution spectroscopy of Ha emission line profiles for 
young low-mass objects, identify likely accretion candidates, and estimate mass accretion 
rates through profile modeling. This work extends measurements of accretion to the lowest 
masses yet examined, approaching the planetary mass regime {M/Mq ~ 0.015), and sig- 
nificantly increases the number of substellar objects with accretion measurements. When 
these results are combined with our previous studies of accretion at higher masses, we are 
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able to trace the change in accretion rates from stars to brown dwarfs across two orders of 
magnitude in mass. 

2. Observations 

For this study, we considered spectroscopically confirmed low-mass members of the 
nearest star forming regions that are optically accessible (/ < 20). Preference was given 
to objects with strong Ha emission. The resulting sample consists of 39 members of the 
Taurus and Chamaeleon I star-forming regions taken from Briceiio et al. (2002), Luhman et 
al. (2003a), Luhman (2004a,b,c), and Luhman (in preparation), which are listed in Table 1. 
The youth and membership of all of the objects in Taurus and Chamaeleon I have been 
established in these previous studies using low-resolution spectroscopy, determining spectral 
types from various atomic and molecular features and selecting against main sequence dwarfs 
or giants using gravity-sensitive K I and Na I absorption lines that are indicative of pre- 
main-sequence sources. The membership of some of these objects has been independently 
confirmed by the presence of reddening in their spectra and their positions above the main 
sequence for the distance of each region, which indicate that they cannot be field dwarfs 
in the foreground or the background of the cloud, respectively. The detections of accretion 
presented in this work represent additional evidence of youth, since such strong broad Ha 
emission is not seen in field dwarfs. 

Our targets arc placed on the Hertzsprung-Russell (H-R) diagram in Figure 1 using 
temperatures and luminosities estimated in the studies cited above. According to the evolu- 
tionary models of Baraffe et al. (1998) and Chabricr et al. (2000), ~ 25 of these sources are 
brown dwarfs, one of which has a mass of only ~ 0.01 Mq. On the nights of 2003 December 
30 through 2004 January 2, we obtained spectra of these targets with the Magellan Inamori 
Kyocera Echelle (MIKE) on the Magellan H 6.5 m telescope at Las Campanas Observatory. 
The instrument provided full coverage from 3200 to 10000 A at a resolution near 20,000 with 
a 1" sht. The exposure times ranged from 300 to 3600 s. 

The MIKE data were bias-subtracted and trimmed using standard IRAF^ routines. The 
processed frames were then cleaned of cosmic rays with the crutils package. Since there was 
no signal in the blue camera spectra and the first few orders of the red camera spectra, 
we only extracted one dimensional spectra for the 22 reddest orders, providing an effective 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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wavelength coverage of about 5700 to 10000 A. The spectra were extracted with the apall 
routine in the IRAF echelle package, and wavelength cahbrated using ThAr lamp spectra. 

3. Analysis 

3.1. Observed Ha Profiles 

The observed Ha profiles for our entire sample are shown in Figure 2. Profile shapes and 
strengths closely resemble those from previous observations of young low-mass objects. Emis- 
sion is detected in all of the targets; ~ 60% of them exhibit narrow (FWHM< 200 km s^^), 
symmetric profiles, while the remainder display broader, asymmetric profiles, occasionally 
with superimposed blueshifted or central absorption reversals. The former type of emission 
is likely due to chromospheric activity on the substellar surface, while the latter is probably 
produced in accretion flows channeled from a circumstellar disk by the substellar magnetic 
held (MuzeroUe et al. (2000); M03; WB03). At least three objects, 2M J04554757+3028077, 
2M J04141188+2811535, and 2M J11013205-7718249, show blueshifted absorption compo- 
nents superimposed on the accretion emission profiles, indicating accretion-powered mass loss 
analogous to that seen in higher-mass T Tauri stars. One of these, 2M Jl 1013205-7718249, 
shows an extreme P Cygni profile, with no evident emission at blue velocities ^ — 50 km 
s~^, and may be entirely formed in the wind. Such wind-dominated profiles are occasionally 
seen in higher-mass CTTSs with large mass loss rates, such as DR Tau (e.g., MCH). Since 
mass loss typically scales with accretion (Hartigan et al. 1995), it is somewhat surprising to 
see this behavior in substellar objects, given their generally small accretion rates. Substellar 
objects of this type are thus worthy of further scrutiny. 

Based on a comparison between Ha emission line equivalent widths {W\), Ha full- 
widths at 10% of the peak intensity (T^io), and spectral types, WB03 determined an empirical 
boundary between accretion and chromospheric activity of Wx > 40 A and Vio > 270 km 
s"^ for young low-mass objects with spectral types later than M6. In Figure 3, we examine 
the same issue for our sample, including for the first time objects at M/Mq < 0.05. We 
see a significant number of objects later than M6 with Wx > 40 A, but with Vio below the 
nominal 270 km s^^ threshold adopted by WB03. This has also been noted by Jayawardhana 
et al. (2003b), who adopted a lower limit of Vio > 200 km s~^ based on observations of three 
accretors with spectral types later than M6. The lower line widths are hkely a reflection of 
intrinsically smaller accretion flow velocity flelds due to the smaller object masses and radii; 
the maximum velocity from infall is Vinfaii — \^2GM^/R^, which yields an upper limit to 
infall line widths of about 195 km s~^ for the lowest-mass objects in our sample. In order 
to deflne the accretor subsample, we have adopted a 10% width criterion V^io > 180 km s~^. 
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which is at the upper end of the most obvious break between Ha emission characteristics. 
A total of 17 objects meet this criterion. 

We may be missing some objects with broad accretion components that are undetectable 
given the signal-to-noise of our data. Using the accretion models described below, we predict 
a minimum measurable mass accretion rate of ~ lO"-*^^ MQyr"^ for objects with M ~ 
0.05 Mq (this limit will vary somewhat with the mass of the central object because of the 
change in the underlying photospheric continuum). Thus, some of the objects not identified 
as accretors in this work may in fact still be accreting, but below this extremely low rate. 
However, given the significant separation of two populations apparent in the 10% width 
distribution, as opposed to a continuous distribution, we believe the number of missing 
accretors is likely to be small. 

3.2. Model Ha Profiles 

Our goal in this paper is to determine mass accretion rates for our substellar sample. 
There are only two primary indicators currently available for this purpose: measuring veiling 
continuum excess or modeling Ha emission profiles. For objects near or below the substellar 
hmit, the former method has been apphed in just a few cases by WB03 and M03. Unfortu- 
nately, veiling is unmeasurable for the large majority of very low mass objects, given their 
generally tiny accretion rates (M03). Thus, modehng of Ha is the most widely applicable 
primary diagnostic, and we use this method here. 

Previously, we have developed radiative transfer models of permitted emission line pro- 
files for both solar-type CTTSs (MCH) and their low-mass stellar and substellar counterparts 
(M03). These models treat line emission produced by magnetospheric accretion flows, where 
gas accreting through a circumstellar disk is channeled onto the star by the stellar magnetic 
field. The model assumes an ideal dipolar field geometry, and a specified gas temperature 
distribution and mass accretion rate; the gas density distribution is determined directly from 
the geometry and the gas velocity as a function of the mass and radius of the star. We defer 
most of the theoretical details to MCH and M03. 

In brief, for this paper we have adopted the same accretion model, using several values 
of the stellar mass (0.025 - 0.15 Mq) and radius (0.25 - 1 Rq) that are representative of 
our sample (see Table 1). The actual mass range will change somewhat depending on which 
evolutionary tracks are used; however, this should not affect our modeling applications since 
the model gas velocities, and, hence the line widths, are only sensitive to -\/M, and the 
derived masses are unlikely to be off by more than a factor of two. We employ the same 
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temperature constraints as those found to reproduce various emission profiles and line ratios 
for CTTSs with well-determined accretion rates from other observations. As we discussed 
in M03, this assumption should be adequate since the gas is probably heated mechanically 
through processes independent of spectral type. The temperature constraints result in an 
inverse proportionality between the gas temperature and density. Due to the lack of opacity 
broadening in any of the observed Ha profiles, we can rule out M > 10~^ M0yr~^ (M03). 
For lower values of M, the gas temperature can be fairly well- constrained; for T > 10, 000 
K, the gas is already completely ionized, so the line flux will not change appreciably, while 
below this level, there is not enough line optical depth to produce emission at observed 
levels. Further details on the parameter constraints are found in M03; here, we apply the 
same methodology. 

We consider the likely accretors identified in the previous section that have spectral 
types M6 and later. The wind-dominated object 2M Jl 1013205-7718249 lacks emission in 
the blue wing, thus we cannot reliably match an accretion model and do not include it in 
the following analysis. For each of the remaining accretor Ha profiles, we constructed a 
model by adopting the magnetospheric temperature and size constraints discussed in M03, 
using the model substellar mass and radius closest to the empirically determined values, and 
finally varying only the mass accretion rate and inclination angle to find the best match 
to the observed profile shape and flux. A range of accretion rates was explored, lO"-*^^ < 
M < 10~^° M0yr~^. The lower limit corresponds to the level where accretion emission is 
produced at an undetectable level, while values above the upper limit would produce too 
much flux both in the wings (from opacity broadening) and in the line overall compared to 
the data in our sample. The inclination angle was varied from pole-on to edge-on orientations 
of the star/disk system. The comparisons were made by eye; in the cases where the profile 
exhibited central or blueshifted absorption components, only the high-velocity line wings 
unaffected by absorption were considered. 

Our best model matches are shown in Figure 4, with parameters listed in Table 2. 
The accretion rates for these objects and the limits for the remainder of our sample are also 
provided in Table 1. We do not include the results for KPNO 4 as the best model match does 
not agree particularly well in terms of profile shape - the observed profile is quite symmetric, 
with a central absorption reversal strongly suggestive of chromospheric emission. The Vio 
width and EW for this object are marginal for accretion, and could also be explained by 
chromospheric emission broadened by rapid rotation; unfortunately, our spectrum does not 
have enough S/N to measure I'sim. 

As mentioned in M03, there is a systematic uncertainty in our model-derived M values 
of about a factor of 3-5. This is mainly due to uncertainties in the size of the accretion flow. 
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which is not well-constrained for our particular sample. However, observations of infrared 
excesses around brown dwarfs in general infer inner disk hole sizes of roughly 2 — 5 R^,, 
similar to our adopted value (Natta et al. 2002; Liu et al. 2003; Mohanty et al. 2004). 

4. Discussion 

Our results more than triple the number of known substellar accretors (those with 
spectral types later than M6) to 16, and extend measurements of mass accretion rates down 
to the lowest masses yet identified. The new accretion rates for our sample of low-mass 
stars and brown dwarfs, combined with previous estimates for other brown dwarfs and more 
massive stars, arc plotted as a function of mass in Figure 5. The accretion rates for substellar 
objects arc extremely small, mostly < 10^^'^ MqJY^'^, several of which are among the lowest 
values measured to date at M = 4 — 5 x 10^^^ MqJt^^. These data continue a trend 
that we have noted previously for the more massive objects in which mass accretion rates 
depend steeply on mass; here, we find a correlation of M oc M^-^. In addition, the change in 
accretion rate with mass is fairly continuous across six orders of magnitude in M and more 
than two orders of magnitude in mass, which suggests that the same formation mechanism 
is responsible for both the low- and high-mass objects. Folding in previous work mentioned 
above, and our own studies of intermediate mass young stars (Calvet et al. 2004), there is 
now significant evidence that the disk accretion paradigm applies to a range of over two 
orders of magnitude in mass, from 0.02 — 3 Mq. Thus, free-fioating stellar and substellar 
objects in this mass range, comprising all but the most massive O and B stars, likely form 
from the collapse of a molecular cloud core onto an accretion disk and central object. 

There is considerable scatter in M for a given mass all across the mass spectrum, as 
has been noted previously for the higher-mass objects. Most of this scatter probably refiects 
intrinsic differences from source to source, and may be due to age differences, since viscous 
evolution expects a general decrease in the accretion rate with time (Hartmann et al. 1998), 
or differences in initial conditions. We note that the accretion rate measurements for low 
mass members of Ophiuchus from Natta et al. (2004) appear to be larger on average by 
about a factor of 5-10 than those in Taurus, Chamaeleon I, and IC 348. This may be an 
indication of the slightly younger age of Ophiuchus, or possibly different initial conditions 
such as the larger envelope infall rates that characterize protostellar objects in that region. 

Our results also imply that objects down to masses of 0.02 Mq harbor significant, 
ordered magnetic fields. Accretion infall cannot occur unless there is a magnetic field of 
sufficient strength to truncate the inner disk and channel material onto the surface of the 
substellar object. Given the model-derived mass accretion rates, inner disk hole size implied 
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by our adopted inner magnetospheric radius, and empirical mass and radius estimates, we 
predict magnetic field strengths on the order of 100-200 G for the lowest-mass accretors in 
our sample, using the theoretical relations given by Konigl (1991). Significant magnetic fields 
are also likely to exist on the nonaccretors, given their signatures of chromospheric activity 
similar to that seen in M and L field dwarfs (e.g. Mohanty & Basri 2003). 

Finally, given the object masses and typical ages of ~ 1 Myr, the mass accretion rates 
in Figure 5 will have little or no effect on the remaining history of the young stars and brown 
dwarfs we have considered. Presumably, most of the mass is accreted in the earliest phases 
of evolution, cither in the initial collapse of the cloud core or via outbursts initiated by 
disk instabilities. As such, it will be important to eventually identify and study embedded 
protostellar counterparts at the lowest masses. 

We thank the staff at Las Campanas Observatory for their support of these observations. 
K. L. was supported by grant NAG5-11627 from the NASA Long-Term Space Astrophysics 
program. 2MASS is a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science Foundation. 
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Table 1. Data for Observed Sample 
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11 


53 


5 


M6 


25 


2962 





019 





53 





075 


-260.0 


362 


-10 


CHSM 7869 


11 


06 


32 


77 


-76 


26 


21 


1 


M6 




2990 





028 





39 





070 


-72.0 


300 


-10 


KPNO 14 


04 


33 


07 


81 


+26 


16 


06 


6 


M6 




2990 





11 


1 


2 





10 


-14.0 


140 


<-12 


KPNO 3 


04 


26 


29 


39 


+ 26 


24 


13 


8 


M6 




2990 





020 





53 





080 


-144.6 


290 


-10 


ISO 252 


11 


10 


41 


42 


-77 


20 


48 


1 


M6 




2990 





027 





62 





080 


-173.1 


274 


-10 


2M ,711070324-7610565 


11 


07 


03 


24 


-76 


10 


56 


6 


M6 




2990 





013 





43 





070 


-9.3 


169 


<-12 


2M J11080234-7640343 


11 


08 


02 


35 


-76 


40 


34 


4 


M6 




2990 





024 





58 





080 


-5.0 


120 


<-12 


2M J11075993-771B317 


11 


07 


59 


93 


-77 


16 


31 


8 


M5 


75 


3024 





067 





95 





12 


-8.4 


178 


<-12 


2M J11072443-7743489 


11 


07 


24 


44 


-77 


43 


49 





M5 


75 


3024 





042 





75 





11 


-8.8 


114 


<-12 


2M .111173792-7646193 


11 


17 


37 


93 


-76 


46 


19 


4 


M5 


75 


3024 





013 





42 





080 


-11.8 


125 


<-12 


2M .104554801+3028050 


04 


55 


48 


01 


+30 


28 


05 





M5 


6 


3014 





014 





44 





075 


-14.0 


150 


<-12 


2M .110561638-7630530 


10 


56 


16 


38 


-76 


30 


53 





M5 


6 


3014 





031 





64 





10 


-41.4 


283 


-10.8 


2M .110580597-7711501 


10 


58 


05 


98 


-77 


11 


50 


1 


M5 


25 


3091 





018 





47 





10 


-5.0 


137 


<-12 


2M J04554757+3028077 


04 


55 


47 


57 


+30 


28 


07 


7 


M4 


75 


3161 





098 


1 








20 


-24.5 


320 




T23 


11 


06 


59 


07 


-77 


18 


63 


6 


M4 


25 


3138 





12 


1 


2 





18 


-7.0 


247 




T3W 


10 


56 


59 


73 


-77 


24 


39 


9 


M3 


5 


3860 





27 


1 


2 





73 


-9.6 


259 




CHXR49 NB 


11 


11 


54 


00 


-76 


19 


31 


1 


M2 


5 


3430 





34 


1 


7 





44 










T3B 


10 


55 


59 


73 


-77 


24 


39 


9 


MO 


5 


3342 













35 


-8.9 


247 





Note. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 
^2MASS Point Source Catalog. 

^Briceno et al. (2002), Luhman et al. (2003a), Luhman (2004a), Luhman (2004b), Luhman (2004c), and Luhman (in preparation). 
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Tabic 2. Htt Model Parameters for Accretor Siibsample 



ID 


M, (Mo) 


R* [Rq) 


Te// (K) 


in 


logM (Moyr-i) 


KPNO 12 


0.025 


0.25 


2600 


45 


-11.4 


KPNO 6 


0.025 


0.25 


2600 


60 


-11.4 


KPNO 7 


0.025 


0.25 


2600 


75 


-11.4 


2M J04414825+2534304 


0.025 


0.25 


2600 


60 


-11.3 


Cha Ha 1 


0.025 


0.25 


2600 


45 


-11.3 


2M J04390396+2544264 


0.025 


0.25 


2600 


70 


-11.3 


2M J04381486+2611399 


0.05 


0.5 


3000 


60 


-10.8 


Cha Ha 11 


0.05 


0.5 


3000 


30 


-11 


CFHT 4 


0.025 


0.25 


2600 


55 


-11.3 


ISO 217 


0.05 


0.5 


3000 


65 


-10" 


2M J04141188+2811535 


0.05 


0.5 


3000 


55 


-10 


CHSM 7869 


0.15 


1.0 


3000 


85 


-10 


KPNO 3 


0.15 


1.0 


3000 


60 


-10 


ISO 252 


0.15 


1.0 


3000 


45 


-10 


2M J10561638-7630530 


0.05 


0.5 


3000 


60 


-10.8 



Note. — All models calculated with magnetospheric radii Rmag — 2.2—3 Rq and maximum 
temperature Ty^ax — 12, 000 K. 

'^Calculated with modified temperature distribution; see text. 
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Fig. 1. — H-R diagram for the targets of the high-resohition spectroscopy in this work shown 
with the theoretical evolutionary models of Baraffe et al. (1998) {M/Mq > 0.1) and Chabrier 
et al. (2000) {M/Mq < 0.1), where the mass tracks {dotted lines) and isochrones {solid lines) 
are labeled in units of Mq and Myr, respectively. 



-14- 



normalized flux 




Fig. 2. — Observed Ha emission line profiles. The spikes close to -500 and +250 km/s seen 
in some spectra (e.g. J04381486+2611399 and KPNO-2) are sky lines that could not be 



- 15 - 



normalized flux 



: . • • . *.. o o o — — u. • . 





Fig. 2. — Continued. 
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— Continued. 
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Fig. 3. — Equivalent widths of Ha vs. the 10% widths of Ha for our sample of young stars 
with spectral types M6 and later (M/Mq < 0.08). The horizontal dotted line indicates the 
W(Hq;) = 40A nominal limit between weak-lined T Tauri stars and Classical T Tauri stars 
for spectral types M6-7.5, as defined by White & Basri (2003). The vertical dashed line is 
our adopted threshold separating accretors and chromospherically active stars. 
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Fig. 4. — Comparisons of observed (solid) and model (dashed) Ha emission line profiles for 
objects in our sample that have profiles indicative of accretion (Table 2). 
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Fig. 5. — Mass accretion rate as a function of mass. Filled circles are from GuUbring et al. 
(1998), White & Ghez (2001), WB03, M03, Calvet et al. (2004), and Natta et al. (2004), 
and open circles are from this work. 



